arXiv:hep-ph/9507369vl 21 Jul 1995 


NONDEGENERATE FERMIONS IN THE 
BACKGROUND OF THE SPHALERON BARRIER 


Guido Nolte 

Fachbereich Physik, Universitat Oldenburg, Postfach 2503 
D-26111 Oldenburg, Germany 

Jutta Kunz 

Fachbereich Physik, Universitat Oldenburg, Postfach 2503 
D-26111 Oldenburg, Germany 
and 

Instituut voor Theoretische Fysica, Rijksuniversiteit te Utrecht 
NL-3508 TA Utrecht, The Netherlands 

Burkhard Kleihaus 

Fachbereich Physik, Universitat Oldenburg, Postfach 2503 
D-26111 Oldenburg, Germany 

February 1, 2008 


Abstract 

We consider level crossing in the background of the sphaleron barrier for nonde¬ 
generate fermions. The mass splitting within the fermion doublets allows only for an 
axially symmetric ansatz for the fermion fields. In the background of the sphaleron 
we solve the partial differential equations for the fermion functions. We find little 
angular dependence for our choice of ansatz. We therefore propose a good approxi¬ 
mate ansatz with radial functions only. We generalize this approximate ansatz with 
radial functions only to fermions in the background of the sphaleron barrier and 
argue, that it is a good approximation there, too. 
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1 Introduction 


The explanation of the observed baryon asymmetry of the nniverse represents a challenging 
problem. Althongh far from solving this highly complex problem, we know at least what 
featnres a theory must have to allow for an explanation. It is therefore remarkable that 
the standard model fulhlls all three Sakharov-conditions to generate the observed baryon 
asymmetry: C and CP violation, a hrst order phase transition and non-conservation of 
baryon number [|l|. 

Here we are concerned with the violation of baryon number (or more generally fermion 
number) in the standard model. It was discovered by’t Hooft 0 as a consequence of the 
Adler-Bell-Jackiw anomaly present in chiral gauge theories. In particular’t Hooft studied 
the fermion number violation induced by vacuum to vacuum tunneling processes described 
by instantons, resulting in extremely small tunneling rates. 

In Weinberg-Salam theory topologically distinct vacua are separated by hnite energy 
barriers. The height of the barriers is given by the energy of the sphaleron, an unstable 
solution of the static field equations Q. Thus the sphaleron determines the minimal 
energy needed for a classically allowed vacuum to vacuum transition. The probability for a 
transition is expected to be enhanced signihcantly, if enough energy is put into the system 
under consideration, either in suitable (future) accelerators or at high temperatures in the 
early universe [5-10]. 

While the barrier is traversed baryon number violation may be seen explicitly by 
analyzing the corresponding Dirac equation in the bosonic background helds. The lowest 
positive energy continuum state becomes continuously deformed along the barrier until 
it reaches the negative energy continuum, passing zero energy precisely at the top of the 
energy barrier, at the sphaleron [11-15]. Investigating the whole spectrum of the Dirac 
equation shows, that along the barrier in fact all levels become continuously deformed into 
the next lower levels, resulting finally in an identical spectrum, where only the number of 
occupied levels above the Dirac sea has decreased by one 
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These calculations [11-16] are based on the approximation, that the fermion doublets 
are degenerate in mass (and that the Weinberg angle may be set to zero BUI). allowing 
for a spherically symmetric ansatz for the fermion wave function. For the physical situa¬ 
tion of highly nondegenerate fermion masses (at least for the heavy flavours) an analogous 
calculation is far more involved, since the spherically symmetric ansatz fails and the equa¬ 
tions of motion cannot be reduced to ordinary differential equations. (This is in contrast 
to the case of instantons 0.) 

Here we consider an axially symmetric ansatz for the fermion helds in the background 
of the sphaleron barrier. The ansatz is chosen in such a way, that it is “almost spherically 
symmetric”, in the sense that the functions involved have little angular dependence. Due 
to the symmetry of the sphaleron the ansatz simplihes considerably in the background 
held of the sphaleron. In this case we solve numerically the full set of partial diherential 
equations for the fermion functions. We then consider a set of approximate ordinary dif¬ 
ferential equations for the fermion functions, hnding almost identical solutions. Because 
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of the numerical complexity involved in solving the full set of partial differential equa¬ 
tions in the background of the sphaleron barrier, we consider in this general case only an 
approximate set of ordinary differential equations for radial fermion functions. We argue 
that these equations represent a good approximation as well. 

In section 2 we briefly review the Weinberg-Salam Lagrangian (for vanishing mixing 
angle) for nondegenerate fermion doublets. In section 3 we present our axially symmetric 
ansatz for the fermions, constructed as a generalization of the usual spherically symmetric 
ansatz. In section 4 we consider fermions in the background of the sphaleron. We derive 
the equations of motion, present the solutions of the full set of partial differential equations, 
and compare with the solutions of the set of approximate ordinary differential equations. 
In section 5 we consider fermions in the background of the sphaleron barrier. We present 
our conclusions in section 6. 


2 Weinberg-Salam Lagrangian 

We start with the bosonic sector of the Weinberg-Salam theory in the limit of vanishing 
Weinberg angle, where the electromagnetic held decouples and can be set to zero, 

Cy, = - A(4>t$ - (1) 

with the held strength tensor 

C = - d.V; + , (2) 

and the covariant derivative 

\igT-V; . (3) 

The SU(2 )l gauge symmetry is spontaneously broken due to the non-vanishing vacuum 
expectation value v of the Higgs held 

( 4 ) 

leading to the boson masses 

Mw = Mz = ^gv , Mh = vV^ . (5) 

We employ the values Mw = 80 GeV, g = 0.65. 

For vanishing mixing angle, considering only one fermion doublet, the fermion La¬ 
grangian reads 


Ff = qLiYDf^qL + qui^d^^qn 

- f^'^^qL^UR + UR^'^qL) - f^‘^\dR^'^qL + qi^dR) , ( 6 ) 
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where qi denotes the lefthanded doublet {uL,dL), while qr abbreviates the righthanded 
singlets {uRjdn), with $ = fr 2 ‘h*. The fermion masses are given by 




(7) 


The fermion equations read in dimensionless coordinates (chosen in units of Mw) 


and 


d d 1 

{i— + *0-*^ + + ^mMTz)qR = 0 


d d 

{i— - ia'-—)qR - {mM^ + AmTzM^)qL = 0 , 


where M is the Higgsheld matrix dehned by 




( 8 ) 

(9) 


( 10 ) 


and m and Am are the average fermion mass and half the mass difference (in units of 

Mw) 

m = {Mu +Md)/{2Mw) : (11) 

Am = {Mu — Md)/{2My/) . (12) 


3 Ansatz 


For the gauge and Higgs fields along the sphaleron barrier we take the usual spherically 
symmetric ansatz in the temporal gauge 


* I 

ya 

*^0 

$ 


l-fA{r)^ . , fB{r)^^ , fc{r)^^ 

T (Aa f if a) T f if a i 

gr gr gr 

0 , 

-^[H(r)+iT- fK(r)) ( J ■ 


( 13 ) 

( 14 ) 

( 15 ) 


Due to a residual gauge degree of freedom we are free to choose the gauge fc = 0. 

To construct an appropriate ansatz for nondegenerate fermions we begin by recall¬ 
ing the spherically symmetric ansatz for degenerate fermions with Am = 0 [11-16,20], 
containing four radial functions, 

qL{r, t) = e“*‘^*M,|.[GL(r) id ■ rFi(r)]xh , (16) 

gnir, t) = -h id ■ fFR{r)]xh , (17) 
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where the normalized hedgehog spinor Xh satishes the spin-isospin relation 

^Xh + rxh = 0. (18) 

The generalized axially symmetric ansatz contains the spherically symmetric ansatz, where 
the four functions Gl, Fl, Gr and Fr now depend on the variables r and 9. Because of 
the presence of the r^-terms in the held equations (||)-(H) for Am ^ 0, we need to ‘double’ 
the ansatz by adding terms of the same structure, but with Xh replaced by TzXt, involving 
the four new (r and 6*-dependent) functions AGr, AFr, AGr and AFr. The ansatz now 
contains eight functions, which are in general complex and ^-dependent, caused by various 
occurrences of the nonvanishing anticommutator [r ■ f,Tz]+ = 2cos6' in the equations of 
motion. Considering the 6'-dependence of the functions, the real part is even in cos 6 while 
the imaginary part is odd. This then suggests the following parametrization of the general 
axially symmetric ansatz, involving 16 real functions of the variables r and p = cos^ 6, 

qR{f,t) = e~"‘^^M^(^[Gl{r,p) + i cos{e)Gl{r, p) + ia ■ f{Fl{r,p) + icos{9)Fl{r,p))] 
+Tz[AGl{r,p) +icos{e)AGl{r,p) 

+ia ■ r{AFl{r,p) + icos{9)AFl{r,p))])xh , (19) 


qR{r, t) = ([G)j(r, p) + i cos(0)G'^(r, p) + ia ■ r(F^(r, p) + i cos(0)F^(r, p))] 

+Tz[AG]i{r,p)+icos{9)AGl{r,p) 

+ia ■ f{AF^{r,p) + icos{e)AF^{r,p))]^Xh ■ (20) 


The choice of ansatz (p!^ -(pOD is not unique. We have also considered alternative 
parametrizations of the axially symmetric fermion ansatz. These involve different fermion 
functions, uniquely related to the above fermion functions. The crucial advantage of the 
ansatz (p!^ - (^0D lies in the observation, that its fermion functions have only a very weak 
angular dependence in the background held of the sphaleron, as shown below. This is in 
contrast to the alternative parametrizations considered. 


4 Sphaleron 

We hrst consider fermions in the background of the sphaleron. Since the background 
held barrier is symmetric about the sphaleron, the fermion eigenvalue is precisely zero at 
the sphaleron [11-16], also for nondegenerate fermion masses. As for degenerate fermion 
masses, the fermion ansatz (©-(Ei then simplihes signihcantly in the background held 
of the sphaleron. This is due to the parity rehection symmetry of the sphaleron, for which 
the functions Jr and H vanish, resulting in the decoupling of eight of the 16 functions. 
These functions, Fl, Gl, AFl, AGl and n, Gl, AFl, can therefore consistently 
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be set to zero. After dropping the number index on the remaining eight functions the set 
of partial differential equations in the variables r and p reads 


0 = -GR’ + -p-^GR + -{l + 2p^)AFR-mKGL 
r op r op 

+AmKAGL — 2pKmAFL , ( 21 ) 

0 = -AGr' + \l + 2p^)FR + -p^AGR + mKAGL 
r op r op 

—AmKGL — 2pKAmAFi , ( 22 ) 

0 = + 1(3 + 2p^)Fr + --^AGr - AmKAFL + mK{FL + 2 AGl) , (23) 

r op r op 

0 = AFR' + -{3 + 2p-^)AFR + --^GR-mKAFL + AmK{FL + 2AGL), ( 24 ) 
r Op r op 


0 = Gl'-- p-^GL-hl + 2p-^)AFL + mKGR + AmKAGR 
r op r op 

+2pK{mAFR + AtuFr) + ^^(G^ + pAFl) , (25) 

r 

0 = AGl --P-^AGl--{I+ 2p-^)Fl - K{mAGR +AitiGr) , (26) 

r op r op 

0 = -Fl'--{ 3 + 2p-^)FL---^AGL-mKFR-AmKAFR 
r Op r Op 

-2K{mAGR + AttiGr) + ^-^{Fr + AGr) , (27) 

r 

0 = -AFr' --{3 + 2p-^)AFr---^Gr + AmKFR + mKAFR . (28) 

r Op r Op 


Inspection of the equations shows, that only three equations, eqs. (PD,(^^ and (^5]), 
contain p-dependent terms, when the terms involving the partial derivative with respect 
to p, present in all eight equations, are not considered. In fact only three functions occur 
with a prefactor p. These are Fr, AFr and AFr. If these three functions are small, then 
the ansatz is approximately spherically symmetric in the sense, that all functions have 
little angular dependence. In the following we show, that this is indeed the case. 

Let us denote the three functions Fr, AFr and AFr as b, as ‘bad’ functions, and the 
other hve functions as p, as ‘good’ functions. First we note, that we could set all three 
bad functions b consistently equal to zero, if the source term 


s = -K{Fr + 2AGi 


(29) 


for the bad functions Fr and AFr in eqs. (^) and (^4]) did vanish. Then the hve good 
functions g were pure radial functions. Let us therefore inspect this source term more 
closely and split it into two terms, s = ai — 02 , with 


ai = -KFr , 


(30) 
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and 


0,2 — ‘2KL • 


(31) 


If oi = 02 , the source term vanishes. We now argue that oi and 02 are approximately 
equal. Setting the bad functions Fr, AFr and AFr equal to zero, and neglecting terms 
with prefactors for large r eqs. and (^71) reduce to 


AGr' = K{mAGR + AmGn) , 


and 

Fr' = -2K{mAGR + AtuGr) . 

With the proper boundary conditions at infinity we thus find for large r for the solu¬ 
tions the desired behaviour, Fr = —2AGr, i.e. the source term vanishes there. On the 
other hand, for small r the source term vanishes, since the function K vanishes. In the 
intermediate region the size of the source term needs numerical analysis. 

We have solved the set of partial differential equations in the background of the 
sphaleron numerically for various values of the average mass m and the mass difference 
Am. Let us consider a typical numerical result. In Fig. 1 we show the ‘good’ lefthanded 
functions, Gr, AGr and Fr, with normalization Gr(0) = 1, for three values of the angle 6 
{6 = 0,7r/4 and 7r/2) for the mass parameters m = 0.5 and Am = 0.25. The ^-dependence 
of the functions is too small to be seen in the figure, being on the order of 10“^. The 
corresponding bad lefthanded function AFr is very small, indeed. For the case considered 
it is less then 5 ■ 10“^, i.e. two orders of magnitude smaller than the good functions, with 
almost no ^-dependence at all. 

These results suggest to approximate all functions by radial functions. We have 
therefore obtained a new set of ordinary differential equations by integrating out the 
^-dependence in the energy density, before variation with respect to the fermion func¬ 
tions. The resulting equations then differ only in prefactors for the three bad functions, 
apart from the absence of the partial derivatives with respect to p. In block form the 
approximate set of differential equations reads 


9' 

b' 


A B 
G D 



(32) 


where A is a 5 by 5 matrix, F is a 5 by 3 matrix etc. The vector Gg represents the 
source terms of the good functions g for the bad functions b. (It is identical in both sets of 
equations.) These source terms are ms, Ams and zero for Fr, AFr and AFr, respectively, 
with the source s defined in eq. (|2D|). 

Solving the approximate set of ordinary differential equations leads to results almost 
identical to those of the full partial differential equations. This is demonstrated in Fig. 1, 
where also the approximate good lefthanded functions Gr, AGr and Fr, with normaliza¬ 
tion Gr(0) = 1, are shown. The difference of the approximate functions and the exact 
functions is too small to be seen in the hgure, being on the order of 10“^. The bad 
lefthanded function AFr is less then 10“^. 
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Thus the exact calculation and the radial approximation result in almost identical 
results, and the bad functions are very small, indeed. We are therefore free to present in 
the following only results obtained with the approximate calculation. In Figs. 2-4 we show 
the same good lefthanded functions, Gl, and as in Fig. 1 for the same value of 
the average mass m = 0.5, but for three different values of the mass difference. Am = 0.25, 
0.5, and 0.75. Fig. 5 is the corresponding figure for the good righthanded function Gr. 
The functions Gi and Gr are the only functions which do not vanish in the limit Am = 0. 
All other functions, which vanish for Am = 0, are approximately proportional to Am as 
seen in Figs. 3 and 4. Finally in Fig. 6 we demonstrate the approximate cancellation of 
the source terms oi and 02 , responsible for the fact that the bad functions are very small. 


5 Sphaleron Barrier 


Let us now consider nondegenerate fermions in the background of the sphaleron barrier. 
Along the barrier we expect a smooth transition of one fermion level from the positive 
continuum to the negative continuum. In the case of degenerate fermion masses, all 
fermion levels change along the barrier to the respective next lower level |^, thus only 
one level crosses zero, and the spectrum exhibits no crossing of any two levels. Expecting 
the same qualitative behaviour of the spectrum in the case of nondegenerate masses, the 
lowest free fermion level, corresponding to the lower mass fermion of the doublet, then 
should cross zero. 


In the general background of the sphaleron barrier the full ansatz, eqs. ([T^)-(p0[), is 
needed. The background helds along the barrier may be taken from the extremal path 
calculations or, as done here, from the gradient approach The set of partial 


equations for the 16 real fermionic functions of the variables r and p reads 


0 = ^Fl-G^^ + ^p^Gi+kl + 2p^)AFl-m(KGl + HFl) 

-Am{-KAGl + HAFl) - 2pKmAFl , (33) 

0 ^ a,Fl- G?,' + i(l + 2pT)G5, - ~AFk - Am(-KAGl + HAFl) 

-m{KGl + HFl - 2KAFI) , (34) 

0 = ojAFl, - AG!,' + 1(1 + 2pA)F2 ^ Ip^AGl - m(-KAGl + HAFl) 

-Am{KGl + HFl) “ 2pKAmAFl , (35) 

0 = u^AFl - AGl' + i(l + 2pAaG% - —Fj, - m(-KAGl + HAFl) 

r op r op 

-Am{KGl + HFl “ ^KAFl) , (36) 

= u^Gl + fH + i(2 + 2pA)Fl - 1(1 + 2p2-)AGl - m(HGl - KFl) 

-Am{HAGl + KAFl) - 2pKmAGl , (37) 


0 











0 - uG\ + + -(3 + 2p—)F^ + AG}j - AmiHAGl + KAFI) 

-m{HGl - KFl - 2KAG\) , (38) 

0 = uAG], + AF}; + -^{2 + 2p^)AF},--^{l + 2p^)Gl-Am{HGl-KFl) 

-m{HAGl + KAFl) - 2pKAmAGl , (39) 

0 = uAG\ +AF'^ +-{?,+ 2p-^AF‘^ +-—G]^-m{HAG\ +KAFl) 

-Am{HGl - KFl - ^KAGl) , (40) 

0 = uFl + Gl'-^p^Gi-\l + 2p^)AFl + m{KG],-HFl)+Am{KAG],-HAFl) 

+2pK{mAFl + AmFl) + ^^{Gl + pAF|) + y (F^ - pAGl) , (41) 

0 = uFl + G]; -\{l + 2p^)Gl+‘^-^^AFl + m{KGl-HFl) + Am{KAGl- HAFl) 

-2K{mAFl + AmFl) + ^^(Gi - AF^) + ^(F| + AG^) , (42) 

0 = o;AFl + AGi'--(l + 2p|-)F|--p|-AGi-m(FAG)j + FAF^) 

r op r op 

-Am(A'G!, + HFi) , (43) 

0 = + AGi'- i(l + 2p|-)AGi +-|-Fl 

r ap r Op 

-m{KAGl + HAFl) - Am{KGl + HFl) , (44) 

"1 1 

0 = c^Gi - Fi' - -(2 + 2 p^)f; + -(1 + 2p^)AG| - + AFA) - Am(ff AG), 

+AAFi) + 2pK(mAGl + AmGj) + ^4d(F2 - pAGl) -KgI + pAFl) , (45) 

0 = ujGI- Fi -)(3 + 2pA)Fl-')-2^AG\-m(HGl + KFl)- Am(HAGl+ KAF^a) 

—2A'(mAG)f + AjnGjf) H- ~~(^l + + ^^A^) , (46) 

0 = ^AGl - AFl' - i(2 + 2 p|-)Af; + 1(1 + 2p|-)Gi - Am(FG;, - KF'g) 

r Op r Op 

-m(HAG\, - KAFp , (47) 

0 = uAGl - AFf' - 1(3 + 2p|-)AF| - —Gl - Am(HGl - KF^) 

r Op r Op 

-m{HAGl - KAFl) ■ (48) 

These equations are analogous in structure to the equations in the sphaleron back¬ 
ground, with all relevant features ‘doubled’. Now six equations contain p-dependent 
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terms (apart from the terms containing partial derivatives with respect to p). These 


are eqs. 


^), (^), (|37D , (|3^ , (1^) and (^5[) . And six functions occur with a prefactor p, 
these are Fl AF|, AG^ and AG|, the six ‘bad’ functions, b. The other ten 

functions are the ‘good’ functions, g. Again, if the bad functions are small, all functions 
have little angular dependence, and an approximation with radial functions only will be 
good. 

Let us therefore inspect the two source terms for the bad functions. 


Si = HGl - KFl - 2KAGI , 


(49) 


S2 = KGl + HFl - 2KAFI , 


(50) 


occurring in eqs. (|38[), (H). and in 
terms according to si = oi — 02 , with 


respectively, and split these two source 


ai = HGl - KFl , (51) 

02 = 2K/\Gl , (52) 

and S 2 = — & 2 , with 

bi = KGl + HFl , (53) 

62 = 2KKFI . (54) 

If both source terms are small, then the bad functions are small, and consequently the 
angular dependence of all 16 fermion functions is small. 

Due to its great complexity, we have not yet attempted to solve the full set of 16 coupled 
partial differential equations numerically. Instead we have from the beginning resorted 
to the study of the approximate set of 16 ordinary differential equations, obtained by 
integrating out the angular dependence in the energy density. But even this approximate 
set of 16 ordinary differential equations has resisted a numerical solution along the full 
sphaleron barrier. Only by setting two of the 16 radial functions explicitly to zero, namely 
the supposedly small bad functions AG| and AG|., we have succeeded in constructing 
the fermion solution along the sphaleron barrier. (Note, that AG| has no source term.) 

Without the solution of the partial diffential equations to compare with, the quality 
of the approximate solution is not known along the full barrier, away from the sphaleron. 
At the sphaleron the approximation is excellent, and it should remain good close to the 
sphaleron. Away from the sphaleron, however, we can at least make a consistency check 
for the radial approximation used, by inspecting the source terms si and S 2 in this ap¬ 
proximation. Numerical analysis shows, that the source terms are indeed small. In Fig. 7 
we show as a typical example along the barrier the source terms bi and 62 for the Chern- 
Simons number Ncs = 0.4 and the mass parameters m = 0.5 and Am = 0.25. While the 
cancellation of the source terms ai and 02 remains as good along the barrier as it is at 
the sphaleron (shown in Fig. 6 ), the cancellation of the additional source terms bi and &2 
is even much better. This indicates, that the bad functions are indeed small compared 
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to the good functions. The radial approximation therefore should be good along the full 
sphaleron barrier. 

Let us then discuss the level crossing along the sphaleron barrier, as obtained with the 
approximate radial set of equations. In Fig. 8 we present the fermion eigenvalue along the 
barrier for an average mass of m = 2 and for several values of the mass difference, Am = 
0.5, 1.0 and 1.5. The eigenvalue starts from the positive continuum at the lower mass 
(1.5, 1.0 and 0.5, respectively), and reaches the negative continuum at the corresponding 
negative value. The bigger the mass splitting, i.e. the smaller the lower mass, the later the 
fermion level leaves the continuum to become bound, analogous to the case of degenerate 
fermion masses [14,15,20]. 

For degenerate fermion masses the fermion wavefunction is determined by the hedgehog 
spinor Xh, giving both isospin components of the fermion doublet an equal amplitude along 
the sphaleron barrier. For nondegenerate fermion masses this is no longer the case. Let 
us dehne the up-part of the fermion wavefunction along the barrier as 

< PT,PT > 

< T, T > 

where P projects out the upper isospin component. (Note, that this dehnition of the 
up-part is not gauge invariant.) For degenerate fermions the up-part is everywhere one 
half. For nondegenerate fermions the up-part along the barrier depends on the size of the 
mass splitting, as shown in Fig. 9 (for the mass parameters employed also in Fig. 8). The 
up-part dominates slightly in the vicinity of the sphaleron and clearly disappears when 
the vacua are reached. Remarkably, the point where the down-part equals the up-part 
only depends on m and not on Am. 

We hnally address the question, how to best approximate a fermion solution in the 
physical situation of nondegenerate fermion masses by a far simpler solution, obtained 
in the approximation of degenerate fermion masses, in the vicinity of the sphaleron. In 
the physical case the mass parameters m and Am determine the nondegenerate masses, 
m -|- Am and m — Am. Close to the continuum clearly the lower mass, m — Am, is the 
relevant fermion mass. In the vicinity of the sphaleron, however, it is the average mass 
m, which matters. In fact, the average mass m of the nondegenerate case mostly leads 
to an excellent approximation for the fermion eigenvalue in the vicinity of the sphaleron, 
when employed in the far simpler calculations with degenerate fermion mass. This is 
demonstrated in Fig. 10, where we compare the nondegenerate case m = 2, Am = 1 with 
the degenerate cases m = 1, Am = 0 and m = 2, Am = 0. Having the same average 
mass, the fermion eigenvalues in the nondegenerate case, and in the second degenerate 
case, agree very well in the vicinity of the sphaleron. In Fig. 11 we present the slope of the 
fermion eigenvalue at the sphaleron as a function of the mass difference, for three values 
of the average mass, m = 0.5, 1 and 2. We observe, that the slope is fairly independent 
of the mass difference Am for not too large values of the average mass m. 
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6 Conclusion 


We have considered level crossing in the background held of the sphaleron barrier for 
fermion doublets with nondegenerate masses. The mass splitting necessitates a gener¬ 
alized ansatz for the fermions, possessing only axial symmetry. We have proposed a 
particular parametrization of the axially symmetric ansatz, containing 16 real functions 
of the two variables r and p = cos^6*. The structure of the ansatz chosen is based on the 
structure of the spherically symmetric ansatz, which represents its simple limit for vanish¬ 
ing mass splitting. This particular parametrization has the great advantage, that it leads 
to fermion functions with little angular dependence in the background of the sphaleron, 
and (supposedly) also along the full sphaleron barrier. 

In the background held of the sphaleron the proposed ansatz simplihes considerably. 
It leads to a set of eight partial differential equations. We have solved these equations nu¬ 
merically, hnding that the resulting fermion functions have very little angular dependence. 
The reason lies in the structure of the equations for this particular choice of ansatz. Only 
three functions occur with an angular dependent prefactor p (apart from partial derivative 
terms), and there is a single source term for these three functions. Since this source term 
is small, these three functions, which introduce explicit angular dependence into the equa¬ 
tions, are small, and consequently all eight functions have only little angular dependence. 

We have then proposed an approximate ansatz with radial functions only. Integrating 
out the angular dependence in the energy density, leads to a new approximate set of 
ordinary differential equations. Solving these numerically, we hnd that the solutions are 
in excellent agreement with those of the full calculation. Thus we have an excellent radial 
approximation for nondegenerate fermion masses at the sphaleron. 

In the general case of fermions in the background of the sphaleron barrier, we have 
found the same structure of the equations as in the sphaleron case, but with all relevant 
features ‘doubled’, since the ansatz no longer simplihes. As yet we have only solved the 
approximate set of ordinary differential equations, obtained by integrating out the angular 
dependence in the energy density (and then setting two of the small functions explicitly 
to zero). Without the solution of the set of partial differential equations to compare with, 
we do not know the quality of the approximation away from the sphaleron. However, we 
have made a consistency check by evaluating the two source terms for the six functions, 
which introduce explicit angular dependence into the equations. Since the source terms 
are small, these functions are small, and consequently all functions have only little angular 
dependence. We therefore argue, that the radial approximation employed should be good 
along the full barrier. 

Considering level crossing along the barrier, we have observed that the fermion mode 
which crosses zero energy at the sphaleron reaches the continua at the lower fermion mass, 
as expected. Finally we have shown, that in the vicinity of the sphaleron the eigenvalue 
for nondegenerate fermions with average mass m and mass difference Am, may be well 
approximated by the eigenvalue obtained with the far simpler calculation, involving only 
degenerate fermions with the average mass m. With respect to the large splitting of the 
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top and bottom quark masses, this suggests to rather use half the top quark mass in 
approximate calculations with degenerate fermions. 
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Sphaleron 



Figure 1: The ‘good’ lefthanded functions, Gl (solid), ^Gl (dotted) and Fl (dashed), 
in the background held of the sphaleron with normalization Gl(0) = 1, in the exact 
calculation for three values of the angle 6 {6 = 0, tt/4 and 7r/2) and in the approximate 
calculation, with the mass parameters m = 0.5 and Am = 0.25. Any of the three visible 
lines consists of four individual lines. 

Sphaleron 



Figure 2: The ‘good’ lefthanded function Gl in the background held of the sphaleron in 
the approximate calculation, for the hxed average mass m = 0.5 and three values of the 
mass diherence Am = 0.25 (solid). Am = 0.50 (dotted). Am = 0.75 (dashed). 
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Figure 3: 

Same as Fig. 2 for 
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Figure 4: 

Same as Fig. 2 for F^. 
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Figure 5: Same as Fig. 2 for Gr. 
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Figure 6: The source term s = ai — a 2 (solid), and its individual parts Oi (dotted) and 
02 (dashed) in the background held of the sphaleron in the approximate calculation, with 
the mass parameters m = 0.5 and Am = 0.25. 
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Figure 7: The source term S 2 = h\ — 62 (solid), and its individual parts hi (dotted) 
and 62 (dashed) in the background held of the sphaleron barrier at the Chern-Simons 
number Ncs = 0.4 in the approximate calculation, with the mass parameters m = 0.5 
and Am = 0.25. 
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Figure 8 : The fermion eigenvalue along the sphaleron barrier in the approximate calcula¬ 
tion, for the hxed average mass m = 2 and three values of the mass difference Am = 0.75 
(solid). Am = 0.50 (dotted). Am = 0.25 (dashed). 
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Figure 9: The up-part of the fermion wavefunction along the sphaleron barrier in the 
approximate calculation, for the fixed average mass m = 2 and three values of the mass 
difference Am = 0.75 (solid), Am = 0.50 (dotted). Am = 0.25 (dashed). 
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Figure 10: The fermion eigenvalue along the sphaleron barrier in the approximate cal¬ 
culation for the average mass m = 2 and the mass difference Am = 1 (solid), compared 
to the fermion eigenvalue for degenerate fermion masses for m = 1 and Am = 0 (dotted), 
and m = 2 and Am = 0 (dashed). 
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Figure 11: The slope of the fermion eigenvalue at the sphaleron in the approximate 
calculation, as a function of the mass difference Am for three values of the average mass, 
m = 0.5 (solid), m = 1 (dotted), m = 2 (dashed). 
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